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CN . Abstract 

ly-v I A general derivation of the muhi-pion correlation function for completely 

chaotic source is given. Its effects on the pion multiplicity distribution, two- 
,,. ^ pion interferometry are studied. A generalized multi-pion correlation function 

■^ I for a partially coherent source is also discussed. 

ix. 

^ ■ I. INTRODUCTION 



C^ 



^ . Hanbury-Brown and Twiss [|l| were the first who apphed the Bose-Einstein correlations 

to measure the size of distant stars. The method was first apphed to particle physics 
by Goldhaber et al. [Q] in 1959. Since then the size of the interaction region has been 
measured by numerous experiments in high energy collisions using many different types 
of particles p|-§l. The basic idea of HBT (GGLP) effect is the following: The intensity 
correlation function of the two identical particles in momentum space is related to the 
Fourier transformation of the space-time distribution of the source. Therefore using the two 
particle correlation function in momentum space, one can extract the radius and lifetime 
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of the source. The source size measurement in high energy heavy-ion colhsions has its 
special meaning: Knowing the size of the source and the stored energy in it one can obtain 
the corresponding energy density reached in the event, which is one of the key quantities 
to investigate the possible transition to quark-gluon plasma (QGP) phase in high energy 
heavy-ion collisions. 

In principle, the extension of such methods to more than two pions is straightforward 
ItHT^. Experimentally, ultrarelativistic hadronic and nuclear collisions provide the envi- 
ronment for creating dozens and in some cases hundreds of pions UlSHTSf. Therefore, it is 



necessary to consider the multi-pion correlations in these processes. The Bosonic nature 
of pion should affect the single and two pion spectrum and distort the two-pion correla- 
tion function. So it is very interesting to analyze the effects of multi-pion correlation on 
these variables 0. Actually the bosonic nature also causes the abundance of pion at small 
momentum. Therefore, people are also interested in analyzing the effects of multi-pion 
Bose- Einstein correlation on the multiplicity distribution |I9|J20|] . 

In 1987, based on Monte-Carlo methods, Zajc analysed the effects of multi-pion cor- 
relation on the phase-space and two-pion interferometry [0. Based on T matrix methods, 
Pratt analysed the effects of multi-pion correlation on the pion multiplicity distribution [ED 



In our paper, we give a more general derivation of the effects of multi-pion correlation on 
two-pion interferometry. The arrangement of this paper is as follows: The derivation of the 
multi-pion correlation function for totally chaotic sources is given in section 2. In section 3, 
the effects of multi-pion correlation on the multiplicity distribution are discussed. In section 
4, the effects of multi-pion correlation on the two-pion correlation function are analysed 
both for N pion events and for events with all possible multiplicities. In section 5, n pion 
correlation function for a partially coherent source is derived. In section 6, as an example, 
our results for a special source distribution are given, together with some discussions. 



II. N-PION CORRELATION FUNCTION FOR A CHAOTIC SOURCE 

The general definition of tlie N-pion correlation function C„(pi, ■ ■ -pn) is 

where Pn{pi ■■■]}„.) is the probability of observing n pions with nionienta{j)i} all in the same 
event. Here we derive a formula for the correlation function which is valid for a totally 
chaotic source. 

A state created by a classical pion source is described by P,pT 



10 >= exp{i dp d'^xj{x)exp{ipx)c^{p})\0 >, (2) 

where c"'"(p) is the pion creation operator. j{x) is the current of the pion, which can be 
expressed as 



here j{x',p) is the probability amplitude of finding a pion with momentum p, emitted by 
the source at x'. A{x') is the probability amplitude of the source distribution. z/(a;') is a 
random phase factor. All sources are uncorrelated in coordinate space when assuming: 

<u*{x')u{x)>=6\x' ~x). (4) 

The coherent state can expand in Fock-Space as 

10 >= E i^iJi^y''c^p)dpdxr^^^^ ^ 1^^^ ^^^ 



Tl' 
n=0,oo ' n=0,oo 



with 



I ^_ {i J j{x)e'P''c-^{p)dpdx)'' ^ 
n\ 

Then we define the n-pion correlation function for n-pion state 

^ /^ ^x <n\c+{pi)- ■■c+{pn)c{pn)---c{pi)\n> 

Cn{Pl, Pn) = pf ,1 , . . , x|, • (7) 



From above definition, the two-pion correlation function for two-pion state can be written 

as 

Cf^ ^ \ — <2|c+(pi)c+(p2)c(p2)c(pi)|2> 



<llc+(pi)c{pi)ll><llc+(p2)c(p2)ll> 



1 I <l|c+(pi)c(p2|l><l|c+(p2)c(pi)|l> 
^ ^ <l|c+(pi)c(pi)|l><l|c+(p2)c(p2)|l> 



(8) 



Noticing that 

c{p)\n >= i / d'^xj {x)exp{ipx)\n — 1 >, (9) 

using eq.(3), we have 

< l\c+{pi)c{p2)\l > = J d'^Xid^X2J*{xi)j{x2)exp{-i{piXi - P2X2)) 

= Jd'x[d^x',fix[,p,)jix',,p2Vix'Mx',) (10) 

A*{x'i)A{x2)exp{—ipix[ +ip2X2). 
Taking the phase average, we get 

< i\c-i.pMmi >= / A;/(.;,p,)7W,P.)|AW)|^e.p(-,to -p,M). (11) 

We define the semiclassical function 

g{x',x,k) = J d'^qj*{x',k + q/2)j{x',k — q/2)exp{iq{x — x')). (12) 

Here k,q are four dimensional momenta. g{x',x,k) can be explained as the probability of 
finding a pion with momentum k at point x due to the source at point x'. Then we have 



j*{x', k + q/2)j{x', k — q/2) = / d'^xg{x', x, k)exp{—iq{x — x')), (13) 

and eq.(lO) becomes 



Here k = ^{pi + P2) , Q = Pi — P2 and g{x,k) is the probability of finding a pion with 
momentum k at point x. It is defined as 

g{x, k) = j d'^x'p{x')g{x', x, k), (15) 

where 

p{x) = A*{x)A{x) (16) 

is the source distribution probabihty. Inserting eq. (11) into eq. (5) gives the correlation 
function 

r(^ -^) = ^ -^ d^xd'^x'gjx, k)g{x', k)exp{iq{x - x')) 

^^''^'^" Jd^xd^x'g{x,pMx',p,) ' ^''^ 

where q = Pi - P2, k = (pi + P2)/2. 

Similarly, the n-pion correlation function can be expressed as 

CniPl,P2, ■ ■ -Pn) = XlXl,<7(l)X2,a(2)---Xn,a(n) (18) 

a 

with the expression of Xij as 

fd*rn(r iPi±PA)piiPi-Pj)^ 

^.^. ^ xiPuPj) = ^ (19) 

^Jd^xd^y{g{x,pi)g{y,pj) 

In the following, we often use the expression of the n-pion momentum probability distri- 
bution Sn which can be expressed as 

Sn = 2j/'l.'^(l)P2,CT{2)---Pn,a(n), (20) 

a 

with 

P^, = P^P^.P^) = - fd'xgix, l^i±M)e^fe-p.)^ (21) 

Here no is the mean multiplicity without B-E correlation, defined as 

no= \ g{x,k)d^xd^k. (22) 



(x{i) denotes the ith element of a permutation of the sequence 1,2,3, ■ ■ -jfi, and the sum 
over (7 denotes the sum over all n! permutations of this sequence. 

Once the source distribution g(x,k) is known, the multi-pion correlation function can be 
calculated from eq.(18). Therefore, C„ can also expressed by S'„ as 

Cn{puP2, ■ ■ -^Pn) = (23) 

Hi Pi,i 

III. EFFECTS OF MULTI-PION CORRELATION ON PION MULTIPLICITY 

DISTRIBUTION 

From the above definition of |n >, we have 



u{n) =< n\n > = ^ < 0|/j*(xi)e '^^'-^^dxi J j{x2)e''-P^''^^dx2dpi 

■ ■ ■ I j*{x2n^l)e'"^"''^"-'dX2n-l I j{x2n)e'P"''^"dX2ndpn\0 >, 

Following a similar derivation as in section two, we get 



(24) 



uj[nj =< n\n > = 

n 



n ,. 

^JS,,l[dp,. (25) 

i=l,n 



Then we have 



k 

<(j)\(j)>=J2<k\k>=J2-^l dp,--- dpkSk. (26) 

k k 

After considering the effect of multi-pion Bose-Einstein correlation, the multiplicity distri- 
bution of pion becomes 

Pin) = <^'^> = ^(^) (27) 

^^ <0|0> EkU^ikY ^'^^ 

For particles without BE Correlation, we have 

Pi,j = ^i,jPi,i (28) 

therefore, 



fSnlldP^ = '^ (29) 

•^ i 

and P{n) reduces to the usual Poisson distribution. From above equation, we can also get 
the mean multiplicity < M > 

<M >=J2n-P{n) (30) 



IV. EFFECTS OF MULTI-PION CORRELATION ON TWO-PION 
CORRELATION FUNCTION 

A. Two-pion correlation function for n-pion events 

For mr events, the two-pion correlation function can be defined as 

where P^(pi,p2) is the modified two-pion inclusive distribution in n pion events. The defi- 
nition of -P"(pi, ■ ■ -.Vi) is 

p„/^ ^^ _ < n\c+{pi) ■ ■ ■ c+{pi)c{pi) ■ ■ ■ c{pi)\n > 

^i ypi) ■ ■ 'Pi) — ~ 1 7 7 • I -I \ • w^j 

< n\n > n ■ ■ ■ [n — t + 1) 

For an event with multiplicity of n, the two-pion inclusive and single-pion inclusive 
distribution can be expressed as 

pn^^ ^^ ^ IUi=3,ndpiSn{pl....Pn) ,^^. 

^ ^' ^ IUi=l,ndpiSn{pi....Pn)' 

and 

pn^^^ IIli=2,ndpiSn{pl....Pn) .^^. 

^ ^ II\i=l,ndpiSn{pi-..-Pn)' 

As n increases, the calculation of the intergration given above becomes more and more 
complex. For the two-pion inclusive distribution in the n pion events, there are only two- 
kinds of terms that we are interested in. One is 
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■ ■ ■ p{Pl,Pk)p{Pk,Pl) ■ ■ ■ p{Pnl,Pn2)p{Pn2,P2)p{p2,Pm)p{Pm,Pn3) " " " p{Pni,Pi)p{PuPl) " -, 

(35) 

the other is 

■ ■■ p{PuPk)p{Pk,Pi) ■ ■ ■ p{Pj,Pi)p{P2,Pi)p{PhPii) ■ ■■p{Pm,P2) ■■■ ■ (36) 

Now we define the function 

Gi{p,q) = nljp{p,pi)dpip{pi,p2)dp2 ■ ■ ■ p{pi^2,Pi^i)dpi^ip{pi-i, q). (37) 

The intergration in above term can be expressed as 

J p{p,Pi)dpip{pi,p2)dp2 ■ ■ ■ dpk-ip{pk-i,q) 

p{q,Pk+i)dpk+ip{pk+i,Pk+2)dpk+2 ■ ■ ■ dpk+i-ip{pk+i-i,p) (38) 

= Gkip,q) ■Gi{q,p) 
and 

J p{p,pi)dpip{pi,p2)dp2 ■ ■ ■ dpi-ip{pi-i,p) 

p{q, Pj+i)dpj+ip{pj+i,Pj+2)dpj+2 ■ ■ ■ dpj+rn-ip{Pj+m~i, q) (39) 

= Gi{p,p)G^{q,q). 
From the expression of S'„ 

Sn = ^ Pl,a{l)P2,a{2)---Pn,a{n) (40) 

the number of terms of the form 

piPl,Pk)piPk,Pj)piPj,Pl) ■ ■ ■ piPm,Pl)piP2,Pil)piPil,Pim)piPim,Pin) " " " piPik, P2)Cn-i 



IS 

(!i^. (41) 

i is the number of p that has not been included in C„_j. Similarly the number of terms of 
the form 



piPl - Pk)p{Pk -Pl)--- p{Pj - P2)p{P2 -Pm)--- p{Pil - Pl)C„ 

is also 

(n-2)\ 



[n — t)\ 



(42) 



Then the two-pion inclusive distribution can be expressed as 



P2{p,q) = —^ TT^~T 5^ [ 51 Gm{P^P)Gi-m{q,q) + Gm{p,q)-Gi-m{q,p)]uj{n~i), 

n{n - 1) uj{n) ,f^„ ^^^_^ 

(43) 
with 

^(^) = 4 / n dpkSn- (44) 

For the single-pion distribution, we are only interested in the terms of the form 
p{Pi,Pk)p{Pk,Pj) ■ ■ ■ p{phPi)Cn~i in the integrand, where i is the number of p that are 
not included in C„_j. The number of this kind terms is 

[n — i)\ 
Then the single-pion distribution is 

PiiP) = -^ E G,{p,pMn - i). (46) 

Now the main question is to extract the expression of uj{n). From the expression of 
eq.(46), we have 

1 

i=l,n 
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uj{n} = - E ^ ■ G{i)uj{n - i) (47) 



with 

C{i) = - fdpG,{p,p). (48) 

t J 

From the above method the two-pion and single pion inclusive distribution can be cal- 
culated for n pion events. 

B. Two-pion correlation function for all events 

For a state |0 > which contains all possible multiplicities, the modified single pion 
distribution can be expressed as 

P'^(-n\ = <'^lc+(p)c(p)|(^> 
-' 1 v^'' «t>\4»<M> 

«t>\^><M> ^ ' 

<4>\4>><M> <M> ' 

where eq.(26) is applied. The two-pion inclusive distribution can be expressed as 

P't'fi^ ^ — <4>\c'^ ip)^^ i^c{g)cip)\4>> 
^2\P-,H) — «t,\<t,><M(M-l)> 

^ Y.U P2iP^9)-nin-l)uj{n) 
~ <4>\4>><M{M-l)> 

(50) 

_ Z]nCi=2,n(I]m=i,i-i ^m (p,p)Gi_„ (g,g)+G„ (p,g)-Gi_™{g,p))w(n-i)] 
~ '' '' <(t>\<t»<M(M-l)> 

_ X:».,^'(P-P)gj'('?-'?)+g'(P''?)-G'i(g'P) 
~ ' <M{M-1)> 

The two-pion correlation function for all multiplicity distribution is 

Pt{q)Pt{p) 

Gi{p,q),C{i) can be calculated either through Monte-Carlo integration or by analytical 
integration. 



cm,i= ^;^)Z,^ . (51) 
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V. N vr CORRELATION FUNCTION FOR A PARTIALLY COHERENT SOURCE 

For a system with one coherent source and many other totally chaotic sources, the state 
is described by 

10 >part= exp{i I dp I d'^x{j{x) + jc{x))exp{ipx)c'^{p))\0 >, (52) 

where c^{p) is the pion creation operator. jc{x) is the current of the pion produced by 
coherent sources, while j{x) is the current produced by totally chaotic sources, which can 
be expressed as eq.(3). The state \(f) >part can be expanded as 

U^ _ V- iliUx) + jUx))e-'P'-'c+{p)dx)'' ^^ ^ _ ^ ,^ ^ 

\(P >part— 2^ j |U >— 2_^ \n >part, W<Jj 

then the definition of the pion correlation function is unchanged. For n-pion state, \n >part, 
the n-pion correlation function can be expressed as 

opart 
Ui=l,n A [Pi) 

with 

Pr\p)=j9{x,p)dx+\UP)\'- (55) 

The S'^"''* can be expressed as 

Sr' = Snipi) + E ^^n,.(p/, Pck) (56) 

k=l,n 

with 

Sn,k{pi^Pc,k) = E Sl{pi,pc,k)- (57) 

;— 1 "! 

'~-^'fe!(n-fc)! 

Here Sn{pi) is similarly defined as eq. (20) but with pj instead of p, the definition of pi is 

PI{P^,PJ) =Jd'xg{x, Pl±Pl)e^(p^-P.^\ (58) 

Sni.PiiPc,k) is similarly defined as Sn{pi) but with k terms pc instead of p/. Certainly, the 
total number of this kind of substitutions is -^{n^z^r I represents one of these substitutions. 
The definition oi pc{pi,Pj) is 

9c{Pi,Pj) =jc{Pi)Jc{Pj)- (59) 
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VI. RESULTS AND DISCUSSIONS 



As an example, we choose the source distribution 



and its non-relativistic form 

where the parameter Rq is the radius, T is the temperature and pg = rriT^T. g{x, ^^) can 
be approximately expressed as 

2 . „2, 2 



P + ^^_ 1 ~k 1 .-^4^x^.^_ 1 -k 1 



— :t- 



g{x,^-—^)= . ^,..,„ e «o- — — e-^2T-^(t) = . .. .,, e % . ,..,„ e % 5(t) 



7ri?g)3/2 (27rm^T)3/2 "^ ^ (Tri?^)^/^ (27rp^)3/2 



(62) 



then we have 



P+^^..(p-.).^^ _ 1 .-^^^^^ 



p(p, g) = y (?(x, £l-:£)eH^-«)-rfx = ^^^^,^3^, 6 ^e ^^0 . (63) 

Define 

Gn{p,q)= p{p,Pi) n dpip{pi,pi+i)dpn^ip{pn-i,q) (64) 

•^ i=l,n-2 

Using eq.(63), we can easily get 

G^{p, q) = a„e-""(p'+^')+^"P^-«^ (65) 



where 



„2 p2 1 p4 

«"+^ = «"-^ = X + A-W (66) 

46„ 4 4pg 166„ 



1 i?^ 

&n = fln + TT + ^ (67) 



4po 4 



^n+1 = 1^ (68) 
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and 



with 



«n+l = a^{^fl\h"' (69) 



R^ 1 1 



Then we have 



^" = ljGn{p,p)dp = «n ,(2a„-,03/2 - (71) 

Now we can easily calculate the multiplicity, single particle and two-particle distribu- 
tion of pions according to the formula given in above sections. In fig.l, we give the pion 
multiplicity distribution. As the BE correlation is included, the probability with high pion 
multiplicity is larger than the one of usual Poissonian distribution. This feature is consistent 
with the nature of Bosons. We also find that as the radius and temperature decreases, the 
effect of BE correlation on the multiplicity distribution becomes larger. 

The parameter no, which is the mean multiplicity without BE correlation, vs. the mean 
multiplicity < M >, is given in fig. 2. When M is not very large, hq increases with increasing 
< M >. However, for very large < M >, with increasing < M > the corresponding no is 



almost a constant. In ref. pO[ to obtain the Poissonian multiplicity distribution for a system 
neglecting symmetrization, a factor -*]- has to be introduced, while in our work, starting 
from the classic multi-pion source the correction of multi-pion BE correlation to multiplicity 
distribution is obtained and it automatically reduces to a Poissonian when BE correlation 
is neglected. 

Including the multipion BE correlation, the single particle momentum distribution for 
events with fixed multiplicity is given in fig.3. It is shown that, as the multiplicity of pion 
increases, the number of pions with low momentum becomes larger. The effect of multi- 
pion correlation on the two-pion correlation function are shown in fig. 4, where we fixed the 
multiplicity of the events. It can be seen clearly that as the multiplicity of the event increases, 
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the two-pion correlation function has a lower chaoticity, though the actual source is totally 
chaotic. The fitted radius also becomes smaller when multiplicity increases as shown in fig. 5. 
This feature can be explained as follows: as the pion multiplicity increases, the emission of 
pions is not independent any more, multi-pion correlation is "Changed into" a factor to 
affect the single-pion emission probability. This gives the source some coherent property 
and makes the pion-source smaller. It can be seen in fig.5 that the coherent parameter 
has some correlations with the radius [^. The lower coherent parameter corresponds to 



the smaller radius. In fig. 6, two-pion correlation function with the mean multiplicity is 
shown. A similar property as in fig.4 can be seen. To obtain the correct behavior of the 
two-pion correlation function at high mean multiplicities it is very important to have correct 
normalization factors in eqs. (49-50), which ensure the correlation function becoming flatter 
and its values approaching to one, when < M > increases. 

As shown in fig. 2, the mean multiplicity increases very fast as no approaches a certain 
value. This seems to show that when < M > becomes very large, the "hardon" picture of 
pion may not be suitable any more, and the system may undergo a kind of phase transition 
into the QGP phase. However, the increase of < M > with increasing n^ should be very 
much controlled by energy conservation which is not included in our formulation. When 
energy conservation is included the behavior of < M > at large hq may change dramatically. 
Therefore, it is worthwhile to further discuss the effects of multi-pion correlation on various 



observables after including energy conservation p3 . 
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Figure Captions 

1. Pion multiplicity distribution with (dot-dashed cure) and without (solid cure) Bose- 
Einstein correlations, the input value of R ,po and no is 3/m, O.lSGeV and 13 respec- 
tively. 

2. no vs. the mean multiplicity M. The dot-dashed line and solid line correspond to 
the case with and without BE correlation. The input value of R and pq is 5fm and 
0.25GeV respectively. 

3. The single particle momentum distribution for different multiplicities, the dashed 
line, dot-dashed line and dotted line corresponds to multiplicity M = 100, 200 and 
400 respectively. The solid line corresponds to the input momentum distribution. The 
input value of R and po is 5fm and 0.25GeV respectively. 

4. Two-pion correlation function for fixed multiplicity events. The solid line, dashed line 
and dot-dashed line correspond to multiplicity M = 100, 200 and 400. The input value 
of R and po is A.95fm and 0.25GeV respectively. 

5. The effective coherent parameter A (dashed line) and radius R (solid line) vs. multi- 
plicity M. The input value of R and po is 4.95fm and 0.25GeV respectively. 

6. Two-pion correlation function for all events. The solid line, dashed line, dot-dashed 
line and dotted line correspond to mean multiplicity < M >= 5, 26, 107 and 180. The 
input value of R and po is 4.95/m and 0.25GeV respectively. 
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